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Abstract 

We examine e + e~ — > DjD* + or D~*Df collisions at 4170 MeV using the CLEO-c detector in 
order to measure the decay constant f D +. We use the Df — > i + v channel, where the £ + designates 
either a u + or a r + . Analyzing both modes simultaneously, we determine B(D+ — > /t/ + ^) = 
(0.657±0.090±0.028)%, B(£>+ -> r+z>) = (7.1 ± 1.4± 0.3)%, and extract / D + = 282±16±7 MeV. 
Combining with our previous determination of B{D + — > u + v), we find that the ratio = 
1.27 ± 0.12 ± 0.03. (All new results here are preliminary.) We compare with current theoretical 
estimates. 



*Submitted to the 33 rd International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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I. INTRODUCTION 



To extract precise information on the size of Cabibbo-Kobayashi-Maskawa matrix ele- 
ments from B — B mixing measurements the ratio of "leptonic decay constants," /$ for Bd 
and B s mesons must be well known Indeed, the recent measurement of B° e — B s mixing 
by CDF [2j that can now be compared to the very well measured B° mixing has pointed 
out the urgent need for precise numbers 0]. The /j have been calculated theoretically. The 
most promising of these calculations are based on lattice-gauge theory that include the light 
quark loops [5( . In order to ensure that these theories can adequately predict fs a / fB d it is 
useful to check the analogous ratio from charm decays f D + / fr,+ . We have previously mea- 
sured fc+ 0, 0| • Here we present the most precise measurement to date of f D + and the ratio 
fD+/fo+- 

The only way in the Standard Model (SM) for a D s meson to decay purely leptonically, 
via annihilation through a virtual W + , is shown in Fig. ^ The decay rate is given by 

F(Dt - tu) = ^f Dt mjM Dt (l - J^-j |Vy 2 , (1) 

where M D + is the Df mass, me is the mass of the charged final state lepton, V cs is a Cabibbo- 
Kobayashi-Maskawa matrix element with a value we take equal to 0.9737 0], and Gp is the 
Fermi coupling constant. 
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FIG. 1: The decay diagram for Dg 



Previous measurements of f D + have been hampered by a lack of statistical precision, 
and relatively large systematic errors 0. One large systematic error source has been the 
lack of knowledge of the absolute branching ratio for — > </>7r + , the mode which most 
measurements have used for normalization |9j. The results we report here will not have this 
problem. 

In this paper we analyze both D s — > and D s — > t + u, t + — > n + V. Both D s decays 
are helicity suppressed because the D s is a spin-0 particle, and the final state consists of 
a naturally left-handed spin- 1/2 neutrino and a naturally right-handed spin- 1/2 anti-lepton 
that have equal energies and opposite momenta. The ratio of decay rates for any two different 
leptons is then fixed by well-known masses. For example, for t + u to the expected ratio 
is 
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T(D+ - r+v) 



Using measured masses j3], this expression yields a value of 9.72 with a negligibly small 
error. After multiplying by £>(t + — > 7r + I7 of 11.06%, the ratio is 1.076 for r + v with respect 
to when r + — > n + T7. 



II. EXPERIMENTAL METHOD 
A. Selection of D s Candidates 

The CLEO-c detector is equipped to measure the momenta and directions of charged 
particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH), detect 
photons and determine their directions and energies [lOj. 

In this study we use 200 pb _1 of data produced in e + e~ collisions using the Cornell 
Electron Storage Ring (CESR) and recorded near 4.170 GeV. Here the cross-section for 
D* + D~+D+D* s ~ is ~1 nb, with D+D~ production being only ~5% of this rate [11]. D 
mesons are also produced mostly as D*D*, with a cross-section of ~5 nb, and in D*D + DD* 
final states with a cross-section of ~2 nb. The DD is a relatively small ~0.2 nb 0|. There 
also appears to be DD production with extra pions. The underlying light quark "continuum" 
background is about 12 nb. The relatively large cross-sections, relatively large branching 
ratios and sufficient luminosities, allow us to fully reconstruct one D s as a "tag," and examine 
the properties of the other. In this paper we designate the tag as a D~ and examine the 
leptonic decays of the D+, though in reality we use both charges. Track selection, particle 
identification, 7, n°, K$ and muon selection cuts are identical to those described in Ref. 

The events we use here occur when e + e~ — > D* + D~ or D+D*~. We will reconstruct 
tags from both final states. The beam constrained mass, tubc-i is formed by using the beam 
energy to construct the mass 0. If we do not detect the photon and reconstruct the m BC 
distribution, we obtain the distribution from Monte Carlo shown in Fig. |21 The narrow peak 
occurs when the reconstructed D s does not come from the D* decay. 

Rather than selecting events based on mBC, we first select an interval that accepts most of 
the events, 2.015 > rriBc > 2.067 GeV, and examine the invariant mass. Distributions from 
data for the 8 modes we use in this analysis are shown in Fig. El Note that the resolution 
in invariant mass is excellent, and the backgrounds not very large, at least in these modes. 
To determine the number of tags we fit the invariant mass distributions to the sum of two 
Gaussians centered at the D~ mass. The r.m.s. resolution (a) is defined as 

o = fiat + (1 - h)a 2: (3) 

where U\ and 02 are the individual widths of the two Gaussians and f\ is the fractional area 
of the first Gaussian. The number of tags in each mode is listed in Table [I] 
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FIG. 2: The beam constrained mass msc from Monte Carlo simulation of e + e 



D+Dt 



(The 



Df — > (j)^ at 4170 MeV. The narrow peak is from the Df and the wider one from the D~ 
distributions are not centered at the Df or D* + masses, because the reconstructed particles are 
assumed to have the energy of the beam.) 



TABLE I: Tagging modes and numbers of signal and background events, within ±2.5<t for all modes, 
except rjp + (±2cr), from two-Gaussian fits to the invariant mass plots. 



Mode 






# 


Background 


K+K-tt- 


h 




8446 ± 160 


6792 


K S K+ 






1852±62 


1021 


T]ir + ; r] — 


, 77 




1101 ±80 


2803 


■q'TT + ; rf - 




4 77 


786±37 


242 


4>p + ; (f) - 


> K+K~, p + - 




1140±59 


1515 


ir + ir~ir + 






2214±156 


15668 


K*+K*°; 


K*+ -» K s ir 




» K-TT+ 1197±81 


2955 


r]p + ; rj -> 


77, p + — > 7T + 




2449±185 


13043 


Sum 






19185 ± 325 


44039 



B. Procedure for Finding Leptonic Decays 

In this analysis we will be selecting events from two processes one where Df —>■ fi + u and 
the other when Df — > r + u, t + — > 7i + V. 1 We first have to select a sample of tag events. We 
require the invariant masses, shown in Fig. El to be within ± 2.5a of the known D~ mass 
(here a is the r.m.s. width). Then we look for an additional photon candidate in the event 



In this paper we use the charge conjugate mode in addition to the specified charge mode. 
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FIG. 3: Invariant mass of D~ candidates (a) K + K~ir + , (b) K$K + , (c) r]ir + , (d) t]'tt + , (e) 4>p + , 
(f) 7r + 7T~7r + , (g) K* + K*°, (h) rjp + , after requiring the total energy to be consistent with the beam 
energy. The curves are fits to two-Gaussian signal functions plus a linear background. 
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that satisfies our shower shape requirement. Regardless of whether or not the photon forms 
a D* with the tag, for real D*D S events, the missing mass squared, MM* 2 recoiling against 
the photon and the D~ tag should peak at the mass. We calculate 

MM* 2 = (E CM -E D - E,) 2 - {-p% - r^) 2 , (4) 

where Eqm is the center of mass energy, En (pd) is the energy (momentum) of the fully 
reconstructed D~ tag, E n (p^) is the energy (momentum) of the additional photon. In 
performing this calculation we use a kinematic fit that constrains the decay products of the 
D~ to the known D s mass and conserves overall momentum and energy. All photons in the 
event are used, except for those that are decay products of the D~ candidate. 

The MM* 2 from the D~ tag sample data is shown in Fig. We fit this distribution to 
determine the number of tag events. This procedure is enhanced by having information on 
the shape of the signal function. One possibility is to use the Monte Carlo simulation for 
this purpose. Our relatively large sample of fully reconstructed D*°D° events allows us to 
examine the signal shape in data when one neutral D is ignored. This sample is shown in 
Fig. |SJ The signal is fit to a Crystal Ball function 0, 14T ]. The a parameter, that represents 



the width of the distribution, is found to be 0.039±0.02 GeV 2 compared with the Monte 
Carlo estimate of 0.310±0.003 GeV 2 . The Monte Carlo does not reproduce well the width 
of the distribution, so we do not use it here. The energy of photons from the D° and D s 
events are somewhat different due to the different masses of the parent hadrons. Thus we 
cannot use the a found here directly. We do, however, get an estimate of the parameters of 
the tail of the distribution. 
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MM* 2 (GeV 2 ) 

FIG. 4: The MM* 2 distribution from events with a photon in addition to the D~ tag. The curve 
is a fit to the Crystal Ball function and a 5th order Chebychev background function. 

Using the fixed tail parameters of the Crystal Ball function, and a 5th order Chebyshev 
polynomial background, we find 12604±423 signal events. After selecting events within an 
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FIG. 5: The MM* 2 distribution from a sample of fully reconstructed D*°D° events where one D° 
is ignored. The curve is a fit to the Crystal Ball function. 

interval 3.978 >MM* 2 > 3.776 GeV 2 , we are left with 11880±399 events. A systematic error 
of 3% is assigned by seeing how the event yields vary when the width of the distribution 
is changed by ±lcr and an additional systematic error of 3% is added in quadrature from 
changing the form of the background function and the fit range. There is also a small 
enhancement of 4% on our ability to find tags in (or r + u, r + — > tt + V) events as compared 
with generic events to which we assign a 10% systematic error or ±0.4%, again added in 
quadrature. 

We next describe the search for Df — > Candidate events are searched for by selecting 
events with only a single extra track with opposite sign of charge to the tag, and also require 
that there not be an extra neutral energy cluster in excess of 300 MeV. Since here we are 
searching for events where there is a single missing neutrino, the missing mass squared, MM 2 , 
evaluated by taking into account the seen D~, and the 7 should peak at zero; the MM 2 
is computed as 



where -Ecm is the center of mass energy, Ejj (pr>) is the energy (momentum) of the fully 
reconstructed D~ tag, E 1 (p^) is the energy (momentum) of the additional photon, and 
(p£) is the energy (momentum) of the candidate muon track. 

We also make use of a set of kinematical constraints and fit each event to two hypotheses 
one of which is that the D~ tag is the daughter of a D*~ and the other that the D* + decays 
into 7-D+, with the Df subsequently decaying into ji + u. The hypothesis with the lowest x 2 



MM 2 = {E cu -E D -E 7 - E^ f - - - p^) 2 , 



(5) 
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is kept. If there is more than one photon candidate in an event we choose only the lowest 
X 2 choice among all the candidates and hypotheses. 
The kinematical constraints are 



Pd s + Pd; = (6) 
-^cm = E Ds + E D * 

E CM , Ml. - Ml a E CM M 2 Dl - M 2 Da 

Ed* = — z 1 777^ or E Da 



2 2£"cm 2 2i?cM 

M D * - M Ds = 143.6 MeV. 

In addition, we constrain the invariant mass of the D~ tag to the known D s mass. This gives 
us a total of 7 constraints. The missing neutrino four-vector needs to be determined, so we 
are left with a three-constraint fit. We preform a standard iterative fit minimizing \ 2 . As 
we do not want to be subject to systematic uncertainties that depending on understanding 
the absolute scale of the errors, we do not make a \ 2 cu t but simply choose the photon and 
the decay sequence in each event with the minimum \ 2 ■ 



III. SIGNAL RECONSTRUCTION 

In this analysis, we consider three separate cases: (i) the track deposits < 300 MeV in the 
calorimeter, characteristic of a non-interacting pion or a muon; (ii) the track deposits > 300 
MeV in the calorimeter, characteristic of an interacting pion; (iii) the track satisfies our 
electron selection criteria defined below. Then we separately study the MM 2 distributions 
for these three cases. The separation between muons and pions is not unique. Case (i) 
contains 99% of the muons but also 60% of the pions, while case (ii) includes 1% of the 
muons and 40% of the pions 0]. Case (iii) does not include any signal but is used later for 
background estimation. 

We exclude events with more than one additional, opposite-sign charged track in addition 
to the tagged D~, or with extra neutral energy. Specifically, we veto events with extra 
charged tracks arising from the event vertex or having a maximum neutral energy cluster, 
consistent with being a photon, of more than 300 MeV. These cuts are highly effective in 
reducing backgrounds. The photon energy cut is especially useful to reject — > 7r + 7r° 
decays, should this mode be significant, and r/7i + decays. 

The track candidates are required to be within the barrel region of the detector | cos#| < 
0.81, where 9 is the angle with respect to the beam. For cases (i) and (ii) we insist that the 
track not be identified as a kaon. For electron identification we require a match between 
the momentum measurement in the tracking system and the energy deposited in the Csl 
calorimeter and the shape of the energy distribution among the crystals is consistent with 
that expected for an electromagnetic shower. 



A. The Expected MM 2 Spectrum 

For the fi + v, final state the MM 2 distribution can be modelled as the sum of two-Gaussians 
centered at zero (see Eq. EJ) A Monte Carlo simulation of the MM 2 for the 07r + subset of 
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K + K 7r + tags is shown in Fig.|B]both before and after the fit. The fit changes the resolution 
er=0.032 GeV 2 to 0.025 GeV 2 , a 22% improvement, and there is no loss of events. 

1000 




.20 



FIG. 6: The MM 2 resolution from Monte Carlo simulation for Df — > fj, + u utilizing a <^>7r + tag and 
a 7 from either D* decay, both before the kinematic fit (a) and after (b). 

We check the resolution using data. The mode Df — > K°K + provides an excellent testing 
ground. We search for events with at least one additional track identified as kaon using the 
RICH detector, in addition to a D~ tag. The MM 2 distribution is shown in Fig. [7| Fitting 
this distribution to a two-Gaussian shape gives a MM 2 resolution of 0.025 GeV 2 in agreement 
with Monte Carlo simulation. 

For the r + u, r + — > 7i + u final state a Monte Carlo simulation of the MM 2 spectra is shown 
in Fig. |HJ The extra missing neutrino results in a smeared distribution. 

B. MM 2 Spectra in Data 

The MM 2 distributions from data are shown in Fig. El Case (i) requires that the candidate 
muon track deposit <300 MeV of energy in the calorimeter, consistent with a non-interacting 
track. Case (ii) requires just the opposite, that the track deposit > 300 MeV and also not be 
consistent with being an electron. Case (iii) requires that the track be consistent with being 
an electron. The overall signal region we consider are below MM 2 of 0.20 GeV 2 . Otherwise 
we admit background from i]tt + and K°ti + final states. There is a clear peak in Fig. EJa), 
due to Df — > yU + z/. Furthermore the region between /x + z/ peak and 0.20 GeV 2 has events 
that we will show are dominantly due to the r + v decay. 

The specific signal regions are defined as follows: for fi + i>, 0.05 >MM 2 > —0.05 GeV 2 , 
corresponding to ±2<r or 95% of the signal; for tv, t + — > 7r + z7, in case (i) 0.20 >MM 2 > 0.05 
GeV 2 and in case (ii) 0.20 >MM 2 > -0.05 GeV 2 . In these regions we find 64, 24 and 12 
events, respectively. 
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FIG. 7: The MM 2 distribution for events with a identified K + track. The kinematic fit has been 
applied. The curve is a fit to the sum of two-Gaussians centered at the square of the K° mass and 
a linear background. 

C. Background Evaluations 

We consider the background arising from two sources: one from real decays and the 
other from the background under the single-tag signal-peaks. For the latter, we obtain the 
background from data. We define side-bands of the invariant mass signals shown in Fig. |2] 
in intervals approximately 4-5<r on the low and high sides of the invariant mass peaks for 
all modes. Thus the amount of data corresponds to approximately twice the number of 
background events under the signal peaks, except for the t]ir~ and r]p~ modes, where the 
signal widths are so wide that we chose narrower side-bands only equaling the data. We 
analyze these events exactly the same as those in the signal peak. 

We list the backgrounds as the number from all modes but t]tt~ and r]p~ first, divided by 
2 and the number for these two modes. For case(i) we find 2/2+1 background in the p + v 
signal region and 3/2+1 background in the r + z/ region. For case (ii) we find 2/2 events. Our 
total background sample summing over all of these cases is 5.5+1.9. 

This entire procedure was checked by doing the same study on a sample of Monte Carlo 
generated at 4170 MeV that includes known charm and continuum production cross-sections. 
The Monte Carlo sample is 7 times the data. We find the number of background events to 
be predicted directly by the simulation to be 28 and the sideband method yields 22. These 
numbers are slightly smaller that what is found in the data but consistent within errors. We 
note that the Monte Carlo is far from perfect in having many estimated branching fractions. 

The background from real Df decays is studied by identifying each possible source of 
background mode by mode. For the p + v final state, the only possible background within the 
signal region is — > 7r + 7r°. This mode has not been studied previously. We show in Fig, mil 
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FIG. 8: The MM 2 distribution for D+ -> t + v, t+ -> ir + V at 4170 MeV. 

the 7T + 7r° invariant mass spectrum. We don't see a signal and set an upper limit < 1.1 x 10~ 3 
at 90% confidence level (C. L.). Recall, that any such events are heavily suppressed by the 
extra photon energy cut of 300 MeV. There are also some Df — > t + u, t + — > n + V events that 
occur in the signal region. We treat these as part of the signal using the Standard Model 
expected ratio of decay rates from Eq. |2] to calculate this contribution. 

For the r + z/, r + — > n + V final state the real backgrounds include, in addition to 
the 7T + 7r° background discussed above, semileptonic decays, possible 7r + 7r 7r° decays and 
other r + decays. Semileptonic decays involving muons are equal to those involving electrons 
shown in Fig. Hfc). Since no electron events appear the signal region, the background from 
muons is also consistent with zero. The 7r + 7r 7r° background is estimated by considering the 
7r + 7r + 7r~ final state whose measured branching ratio is (1.02±0.12)% 0|. This mode has 
larg e contributions from /o(980)7r + and other resonant structures in 7r + 7r~ at higher mass 
[15| . The 7r + 7T 7r° mode will also have these contributions, but the MM 2 opposite to the 
7r + will be at large mass. The only component that can potentially cause background for 
us is the non-resonant part measured by FOCUS as (17±4)%. This background has been 
evaluated by Monte Carlo simulation as well as others from other r + decays, and listed in 
Table El 



IV. CHECKS OF THE METHOD 

We perform an overall check of our procedures by measuring B(Df — > ~K°K + ). For this 
measurement we compute the MM 2 (Eq. EJ) using events with an additional charged track 
but here identified as a kaon. These track candidates have momenta of approximately 1 
GeV/c; here our RICH detector has a pion to kaon fake rate of 1.1% with a kaon detection 
efficiency of 88.5% [16]. For this study, we do not veto events with extra charged tracks or 
photons. 
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FIG. 9: The MM 2 distributions from data using D s tags and one additional opposite-sign charged 
track and no extra energetic showers (see text). For the case when the single track (a) deposits 
< 300 MeV of energy in the calorimeter, case (i). The peak near zero is from Df — > [i^v events, 
(b) Track deposits > 300 MeV in crystal calorimeter but is not consistent with being an electron, 
case (ii). (c) The track is identified as an electron case (iii). 



The MM 2 distribution is shown in Fig. [7| The peak near 0.25 GeV 2 is due to the decay 
mode of interest. We fit this to a linear background from 0.02-0.50 GeV 2 plus a two-Gaussian 
signal function. The fit yields 228±18±8 events. Since r]K + could in principle contribute an 
asymmetric background in this region, we searched for this final state. Not finding any signal, 
we set an upper limit of 2.8xl0 -3 , approximately a factor of five below our measurement. 
In order to compute the branching fraction we must include the efficiency of detecting the 
kaon track 76.2%, including radiation [tH |. and take into account that it is easier to detect 

the photon from the D* decay in ~K°K + events than in the average D S D* event due to the 
track and photon multiplicities, which gives a 4% correction. These rates are estimated by 
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FIG. 10: The invariant 7r + 7r° mass. The curves are fits to a background polynomial and Gaussian 
signal with a width fixed by Monte Carlo simulation. 

TABLE II: Backgrounds in the D+ -> t + v, t+ -> vr+F sample, case (i) for 0.2>MM 2 > 0.05 GeV 2 
and case (ii) for 0.2>MM 2 > -0.05 GeV 2 



Mode 






£>(%) # of events case(i) 


# of events case 


'ii) Sum 


7T+7T 7T 






1.0 


0.1 


0.3 


0.4 


Df - T 






6.2 








T+ -► 






1.6 


0.3 


0.4 


0.7 


r+ -► 






1.1 











r+ -» 






1.1 


0.2 





0.2 


Sum 








0.6 


0.7 


1.3 



using Monte Carlo simulation. We determine 

B(D+ -> K°K + ) = (2.74 ± 0.23 ± 0.13)%, (7) 

where the systematic errors are listed in Table ITLT1 The largest components of the systematic 
errors arise from the number of signal events (4.2%) from the signal function width (3%) 
and the shape of the background function (3%). This method is in good agreement with the 
latest double tag result jl8j |. 

V. LEPTONIC BRANCHING FRACTIONS 

The sum of MM 2 distributions for case (i) and case (ii) normalized to the expectation 
of the sum of — ► fi + u and — > t + u, t + — > n + u is shown in Fig. ^2 The data are 
consistent with our expectation containing mostly signal for MM 2 <0.2 GeV 2 . Recall there 
are 100 total events only 5.5 of which we estimate are background. Above 0.2 GeV 2 other, 
larger backgrounds enter. 
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TABLE III: Systematic errors on determination of the — * K K + branching ratio. 



Error Source 


Size (%) 


Signal shape 


1 


Background shape 


1 


Track finding 


0.7 


PID cut 


1.0 


Number of tags 


4.2 


Total 


4.6 




0.00 



0.25 



0.50 



Missing Mass squared (GeV 2 ) 



FIG. 11: The sum of case (i) and case (ii) MM 2 distributions compared to the predicted shapes for 



7T 1 V. 



The number of fi + u events in the signal region, N^, is related to the number of tags, N tl 
the branching ratios and the estimated backgrounds iVbkgrd as 



Nftv -^bkgrd 



tag 



e'B(Dj - p+v) + e"B(D s 



T ' V, T 



IT' V) 



where e (79.5%) includes the efficiency for reconstructing the single charged track including 
final state radiation (77.8%), the (98.3±0.2)% efficiency of not having another unmatched 
cluster in the event with energy greater than 300 MeV, and for the fact that it is easier to find 
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tags in /i + z/ events than in generic decays by 4%, as determined by Monte Carlo simulation. 
The efficiency labeled e' is a product of the 99% muon efficiency for depositing less than 300 
MeV in the calorimeter and 92.3% acceptance of the MM 2 cut of |MM 2 < 0.05 1 . The quantity 
e" is the fraction of r + z/; r + — > tt + v events contained in the signal window (13.2%) times 
the 60% acceptance for a pion to deposit less than 300 MeV/c in the calorimeter (6.4%). 
The two D + branching ratios in Eq. |H] are related as 

B(D+ -> t + u- t + -> 7T+z/) = R-B(r + -> tt + v)B(D+ -> = 1.076 -> /x + i/) , (9) 

where we take the Standard Model ratio for R as given in Eq. Eland B{t + — > 7r + z/)=11.06% 
0. This allows us to solve Eq. El Since iV^= 64, A^ bkgrd =2, and N tag = 11880 ± 399 ± 499, 
we find 

B(D+ -> = (0.657 ± 0.090 ± 0.028)%. (10) 

We can also sum the /z + z/ and r + v contributions, where we restrict ourselves to the MM 2 
region below 0.20 GeV 2 . Eq. [S] still applies. The number of signal and background events 
changes to 100 and 7.4, respectively. The efficiency e' becomes unity, and e" increases to 
45.2%. Using this method, we find 

B{Df -> = (0.664 ± 0.076 ± 0.028)%. (11) 

The systematic errors on these branching ratios are given in Table IIVI 



TABLE IV: Systematic errors on determination of the — > fi + v branching ratio. 



Error Source 


Size (%) 


Track rinding 


0.7 


Photon veto 


2 


Minimum ionization 


1 


Number of tags 


4.2 


Total 


5.2 



We also can analyze the r + v final state independently. We use different MM 2 regions for 
cases (i) and (ii) defined above. For case (i) we define the signal region to be the interval 
0.20>MM 2 >0.05 GeV 2 , while for case (ii) we define the signal region to be the interval 
0.20>MM 2 >-0.05 GeV 2 . Case (i) includes 98% of the signal, so we must exclude the 
region close to zero MM 2 , while for case (ii) we are specifically selecting pions so the signal 
region can be larger. The upper limit on MM 2 is chosen to avoid background from the tail 
of the K 7r + peak. The fractions of the MM 2 range accepted are 32% and 45% for case (i) 
and (ii), respectively. 

We find 24 signal and 3.1 background events for case (i) and 12 signal and 1.7 backgrounds 
for case (ii). The branching ratio, averaging the two cases is 

B(D+ ^t + u) = (7.1 ±1.4 ±0.3)%. (12) 

Lepton universality in the Standard Model requires that the ratio R from Eq. Elbe equal to 
a value of 9.72. We measure 

* s rU^ r 10 - 8 ^ - 2 - (13) 
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Thus we see no deviation from the predicted value. Current results on D + leptonic decays 
also show no deviations fTg| . The absence of any detected electrons opposite to our tags 
allows us to set an upper limit of 

B(D+ -> e + u) = 3.1 x 1(T 4 (14) 

at 90% confidence level; this is also consistent with Standard Model predictions and lepton 
universality. 

Using our most precise value for B{D^ — > [i^v) from Eq. ^2 that is derived using both 
our /i + z/ and r + z/ samples, and Eq. [TJ we extract 

f Ds = 282 ± 16 ± 7 MeV. (15) 

We combine with our previous result 0] 

f£ = 222.6 ± 16.7±|;1 MeV. (16) 

and find a preliminary value for 

-p- = 1.27 ±0.12 ±0.03, (17) 

JD+ 

where a small part of the systematic cancels in our two measurements. 



VI. CONCLUSIONS 

f 4 

Theoretical models that predict f D + and the ratio y^- are listed in Table IVl Our result 
is higher than most theoretical expectations. We are consistent with Lattice-Gauge theory, 
and most other models, for the ratio of decay constants. 

TABLE V: Theoretical predictions of fjj+ and f D + / fp+ . QL indicates quenched lattice calculations. 



Model 


f D + ( MeV ) 


f D+ (MeV) 


f D +/fD+ 


Lattice (n f =2+l) [20] 


249 ± 3 ± 16 


201 ± 3 ± 17 1.24 ± 0.01 ± 0.07 


QL (Taiwan) [21] 


266 ± 10 ± 18 235 ± 8 ± 14 1.13 ± 0.03 ± 0.05 


QL (UKQCD) [22] 


236 ± 8±{l 


210 ± 10t\l 


1.13 ±0.02jlg;^ 


QL [23J 


231 ± 12^5 


211 ± 14i^ 2 


1.10 ±0.02 


QCD Sum Rules [24] 


205 ± 22 


177 ±21 


1.16 ±0.01 ±0.03 


QCD Sum Rules [25] 


235 ± 24 


203 ± 20 


1.15 ±0.04 


Quark Model [26] 


268 


234 


1.15 


Quark Model [27] 


248±27 


230±25 


1.08±0.01 


Potential Model [28] 


241 


238 


1.01 


Isospin Splittings [29] 




262 ± 29 





By using a theoretical prediction for f D + / f D + we can derive a value for the ratio of CKM 
elements \V cd /V cs \. Taking the value from Ref. of 1.24 ± 0.01 ± 0.07, we find 

\V cd /V cs \ =0.22 ±0.03 , (18) 
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where the theoretical and experimental errors have been added in quadrature. This value is 
consistent with expectations. 

We now compare our preliminary results with previous measurements. The branching 
fractions modes and derived values of f D + are listed in Table IVI1 Our values are shown 
first. We are generally consistent with previous measurements that are also higher than the 
theory. 

TABLE VI: These results compared with previous measurements. Results have been updated for 
new values of the D s lifetime. ALEPH uses both measurements to derive a value for the decay 
constant. 



Exp. 


Mode 


B 


(%) 


f Dt (MeV) 


CLEO-c 




(6.57 ± 0.90 ±0.34)10" 3 




281 ± 19 ± 7 


CLEO-c 


T + V 


(7.1 ± 1.4±0.3)10" 2 




296 ± 29 ± 7 


CLEO-c 


combined 






282 ± 16 ± 7 


CLEO [30] 




(6.2 ± 0.8 ± 1.3 ± 1.6)10~ 3 


3.6±0.9 


273 ± 19 ± 27 ± 33 


BEATRICE [3JJ 




(8.3±2.3±0.6±2.1)10~ 3 


3.6±0.9 


315 ±43 ±12 ±39 


ALEPH [22] 




(6.8 ± 1.1 ± 1.8)10~ 3 


3.6±0.9 


285 ± 19 ± 40 


ALEPH [22] 


T^V 


(5.8 ± 0.8 ± 1.8)10~ 2 






OPAL [34] 


T + V 


(7.0±2.1±2.0)10~ 3 




286 ± 44 ± 41 


L3 [33] 


T + V 


(7.4 ± 2.8 ± 1.6 ± 1.8)10 -3 




302 ± 57 ± 32 ± 37 


BaBar [36] 




(6.5±0.8±0.3±0.9)10~ 3 


4.8±0.5±0.4 


279 ± 17±6± 19 



Most measurements of -D+ — * i + v are normalized with respect to B(D+ — > <f)7i + ). One 
measurement that isn't is that of OPAL, which normalizes to a D s fraction in Z° events that 
is derived from an overall fit to heavy flavor data at LEP j35|. It still, however, relies on 
absolute branching fractions that are hidden by this procedure, and the estimated error on 
the normalization is somewhat smaller than that indicated by the error on available at 
the time of their publication. The L3 measurement is normalized to using a calculation that 
the fraction of D s mesons produced in c quark fragmentation is 0.11±0.02 and that ratio 
of D*/D s production is 0.65±0.10. The ALEPH results use for their /i + z/ results and a 
similar procedure as OPAL for their r + z/ results. We note that the recent BaBar result uses 
a larger B^ than the other results. 

The CLEO-c determination of f D + is the most accurate to date. It also does not rely on 
the independent determination of any normalization mode. 
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